(* Ïåðåäàòî÷íàÿ ôóíêöèÿ ñêîðåêòèðîâàííîé ñèñòåìû äåìïôåðà ðûñêàíèÿ Wy/dí *)
w=(1*7.56 (0.47 +s) (0.47 +s))/((0.47 +s) (41.01 +12.80 s+s2))
Wt=(1*7.56 (0.47 +s) (0.47 +s))/((0.47 +s) (41.01 +12.80 s+s2))
 (7.56 (0.47 +s))/(41.01 +12.8 s+s2)
 (7.56 (0.47 +s))/(41.01 +12.8 s+s2)
 (* Ïåðåéäåì îò s ê êîìïëåêñíîé ôîðìå *)
wk=w/.s®I*omega
 (7.56 (0.47 +ä omega))/(41.01 +12.8 ä omega-omega2)
 f=(Abs[wk])^2
 57.1536 Abs[(0.47 +ä omega)/(41.01 +12.8 ä omega-omega^2)]2
 fvhod=1 (*âõîäíîé ñïåêòð*)
 1
 fvihod=f*fvhod
 57.1536 Abs[(0.47 +ä omega)/(41.01 +12.8 ä omega-omega^2)]2
 (* ÑÏÅÊÒÐ *)

Plot[fvihod,{omega,0,300},PlotRange®Full]
 
 p=NIntegrate[fvihod,{omega,0,Infinity}]
 7.05158
 omega0=0;
shag=0.1;
While[shag³10^-5,{While[(NIntegrate[fvihod,{omega,0,omega0}]-0.95*p)<0,omega0=omega0+shag],omega0=omega0-shag,shag=shag/10}]
omega0=omega0-shag*5
 161.934
 p95=NIntegrate[fvihod,{omega,0,omega0}]
 6.699
 (* ÏÅÐÈÎÄ ÊÂÀÍÒÎÂÀÍÈß *)
Tkvant=0.1*Pi/omega0 (* ãäå 0.1 - ïîïðàâî÷íûé êîýôôèöèåíò*)
 0.00194005
 ht=InverseLaplaceTransform[Wt/s,s,t]
 7.56 (0.0114606 +ã(-6.4-0.223607 ä) t ((-0.00573031+2.07206 ä)-(0.00573031 +2.07206 ä) ã0.447214 ä t))
 Plot[ht,{t,0,10},PlotRange®Full]
 
 hf[t_]=Abs[ht]
 7.56 Abs[0.0114606 +ã(-6.4-0.223607 ä) t ((-0.00573031+2.07206 ä)-(0.00573031 +2.07206 ä) ã0.447214 ä t)]
 (* Âðåìÿ ðåãóëèðîâàíèÿ è ïåðåðåãóëèðîâàíèÿ *)
k=10000;(*êîë-âî ïðîìåæóòêîâ, ãäå íàõîäèì ðåøåíèå*)
dt=10/k;
hust=hf[10000](*óñòàíîâèâøååñÿ çíà÷åíèå*);
hmax=hf[0](*çàäàäèì íà÷àëüíîå ìàêñ. çíà÷åíèå*);
(*×èñëåííîå ðåøåíèå óðàâíåíèé h (t)=0.95 è h (t)=1.05, îáúåäèíåííûõ â îäíî |h (t)-1|-0.05=0 è îïðåäåëåíèå ìàêñèìóìà*)
For[i=1,i<k+1,{If[Sign[Abs[hf[dt*(i-1)]-hust]-0.05]¹Sign[Abs[hf[dt*(i)]-hust]-0.05],Treg=dt*(i-0.5)],If[hmax<hf[dt*i],hmax=hf[dt*i]]};i++]
sigmaper=Round[100*(hmax-hust)/hust,0.01]; (*ïåðåðåãóëèðîâàíèå*)
Print["Treg = ",Treg," ñ"] 
Print["Sigma = ",sigmaper," %"] 
 Treg =  0.7165  ñ
 Sigma =  429.31  %
 (*Êîë-âî ïåðèîäîâ êâàíòîâàíèÿ çà âðåìÿ ðåãóëèðîâàíèÿ*)
Nkvant=Ceiling[Treg/Tkvant]
 370
 Stextrop=3(*Ñòåïåíü ïîëèíîìà ýêñòðàïîëÿöèè*);
Masotklon={}(*Ìàññèâ îòêëîíåíèé*);
hizm=Table[0,{Stextrop}](*Ìàññèâ ïîñëåäíèõ èçìåðåííûõ (ðåàëüíûõ) çíà÷åíèé èñïîëüçóåìûõ ïðè ýêñòðàïîëÿöèè*);
 For[i=Stextrop-1,i<=Nkvant,i++,
  {Ttek=Tkvant*i;
   hTtek=hf[Ttek];
   hizm[[1]]=hf[Ttek-Tkvant*2];
   hizm[[2]]=hf[Ttek-Tkvant*1];
   hizm[[3]]=hf[Ttek-Tkvant*0];
   (*Ýêñòðàïîëÿöèÿ Ìåòîä Àäàìñà*)
   hextrop=hizm[[3]]+Tkvant*(23*hizm[[3]]-16*hizm[[2]]+5*hizm[[1]])/12;
   Masotklon=Append[Masotklon,Abs[hf[Ttek+Tkvant]-hextrop]];
   }];
 (*Ìàòåìàòè÷åñêîå îæèäàíèå îòêëîíåíèÿ èñòèííîãî çíà÷åíèÿ îò ýêñòðàïîëèðîâàííîãî çíà÷åíèÿ*)
Matozh=Mean[Masotklon] 
 0.00234038
 (*Ñðåäíåå êâàäðàòè÷åñêîå îòêëîíåíèå îòêëîíåíèÿ*)
sigma=StandardDeviation[Masotklon] 
 0.00233929
 (* ÏÎÐÎÃ ÑÐÀÁÀÒÛÂÀÍÈß *)
b=3 (*1..5 - äîâåðèòåëüíàÿ ãðàíèöà*);
Porog=Matozh+b*sigma 
 0.00935825